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We study the effects of finite stellar temperatures on the oscillations of superfluid neutron stars. 
The importance of these effects is illustrated with a simple example of a radially pulsating general 
relativistic star. Two main effects are taken into account: (i) temperature dependence of the 
entrainment matrix and (ii) the variation of the size of superfluid region with temperature. Four 
models are considered, which include either one or both of these two effects. Pulsation spectra are 
calculated for these models, and asymptotes for eigenfrequencies at temperatures close to critical 
' temperature of neutron superfluidity, are derived. It is demonstrated that models that allow for 

, the temperature effect (ii) but disregard the effect (i), yield unrealistic results. Eigenfunctions 

for the normal- and superfluid-type pulsations are analyzed. It is shown that superfluid pulsation 
modes practically do not appear at the neutron-star surface and, therefore, can hardly be observed 
by measuring the modulation of the electromagnetic radiation from the star. The e-folding times 
for damping of pulsations due to the shear viscosity and nonequilibrium modified Urea processes 
are calculated and their asymptotes at temperatures close to the neutron critical temperature, are 
obtained. It is demonstrated that superfluid pulsation modes are damped by 1-3 orders of magnitude 
faster than normal modes. 

PACS numbers: 97.60.Jd, 47.75. +f, 97.10.Sj, 47.37.+q 
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^ ■ I. INTRODUCTION 

Q I High-frequency oscillations of the electromagnetic radiation, interpreted as resulting from neutron-star oscillations, 
. have already been observed during giant flares P, Also, detectors are being designed, which, according to prelimi- 
^ ' nary estimates, will be able to detect gravitational radiation from oscillating neutron stars (NSs) To interpret 

the existing and future observations, it is necessary to have a well-developed theory of neutron-star oscillations. Our 
study is devoted to the superfluid sector of this theory. To be more concrete, the goal of the present paper is to 
—i , demonstrate, using radial pulsations as an example, that the effects of finite temperatures strongly affect dynamics 
> ■ of NSs. 

Neutron stars are relativistic objects having supranuclear densities and strongly curving spacetime around them. 
Therefore, accurate calculations of their properties should be made in the frame of the general relativity theory. The 
situation is additionally complicated by the fact that baryons in the neutron-star cores can be in superfluid state 
Thus, several velocity fields can co-exist in matter without dissipation, which makes hydrodynamic equations 
substantially more complicated flol - [T7| . Epstein [31 was the first who considered sound waves in the superfluid matter 
of NSs. Global pulsations of superfluid Newtonian stars were first studied by Lindblom and Mende h I14ll. Since then, 
a growing number of papers devoted to oscillations of superfluid NSs have been published (see, e.g., |l0l - ll3l[l9l - l27| and 
references therein) . To make the problem tractable, most of these authors employed the nonrelativisti c hy drodynamics 
in their calculations. Pulsations of general relativistic superfluid NSs were analyzed only in Refs. [l2|Jl3l [20l - [23| under 
a number of simplifled assumptions. In particular, until recently it was customary to use the zero-temperature limit of 
\^ , superfluid hydrodynamics. However, as it was first shown in Ref. [l^] (hereafter GA06), pulsation spectra of superfluid 
NSs can be very sensitive to variation of temperature. 

The authors of GA06 studied oscillations of superfiuid NSs using the finite-temperature version of relativistic 
hydrodynamics for superfluid mixtures, described in detail in Refs. [ll|, [l^, HO, [3l|. For simplicity, only radial 
pulsations were analyzed and the simplest npe-matter composition of the stellar core was adopted, including neutrons 
(n), protons (p), and electrons (e). Neutrons and protons in the core were allowed to be superfluid, while free neutrons 
in the crust were treated as normal (nonsuperfluid) . Owing to the fact that the relativistic entrainment matrix (analog 



of the superfluid density jl5| for mixtures) is a steep function of temperature [32|, the spectrum obtained in GA06 
was strongly temperature-dependent. The authors of GA06 neglected another temperature effect, the decrease in the 
size of the superfluid reg ion with increasing temperature. It should be noted that this effect has been considered 
in several papers [23l. l28l. [29j : however, the authors of these papers disregarded the temperature dependence of the 
entrainment matrix (that is, they used the zero-temperature limit of superfluid hydrodynamics). 

In the present study we, for the first time, analyze the influence of both of these temperature effects on the pulsation 
spectra. For that we consider the more realistic profiles of the neutron critical temperature Tcn{r) than that employed 
in GA06, where it was assumed that the red-shifted Ten is constant throughout the core. We also compare the 
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eigenfunctions for normal and superfluid pulsation modes, and calculate their characteristic damping times. The 
asymptotes for the eigenfrequencies and for e-folding times at temperatures close to T^n will also be examined. 

The paper is organized as follows. In Sec. II we briefly discuss the superfluid hydrodynamics for mixtures and write 
down a system of equations describing radial oscillations of general relativistic superfluid NSs. In Sec. Ill, we present 
the pulsation spectra calculated for four different neutron-star models that allow for one or both of the temperature 
effects (the temperature dependence of the entrainment matrix and/or the variation of the size of the superfluid 
region with temperature). In Sec. IV, the approximate approach of Refs. (33l.[3^ is used to obtain the asymptotes for 
pulsation eigenfrequencies at temperatures close to Tc„. Sec. V is devoted to analysis of the eigenfunctions of radial 
pulsations. Finally, in Sec. VI we study damping of radial pulsations due to the shear viscosity and nonequilibrium 
modified Urea processes. Section VII presents the summary. 

In the following, unless otherwise stated, the system of units used is one in which the Boltzmann constant ks, the 
Planck constant h, and the speed of light c equal unity: fcs = h = c = 1. 



II. BASIC EQUATIONS 
A. Superfluid hydrodynamics 

Following GA06, we consider the simplest npe-matter composition of the neutron-star cores. Because both protons 
and neutrons can be in the superfluid state, one has to use the relativistic hydrodynamics of superfluid mixtures to 
study the oscillations of NSs (see, e.g., [Ill, HO, The main distinctive feature of this hydrodynamics is the 

presence of several velocity fields in the mixture. In our case, these are the four-velocity of the 'normal' component 
of matter (electrons and Bogoliubov excitations of neutrons and protons) as well as the four-velocities of superfluid 
neutrons v'^^^) ^^^id superfluid protons I'^^p) • Below we briefly discuss the basic equations describing nondissipative 
motion of superfluid mixtures (for more detail, see [iTI. [T2|. Isol. Isij ) . 

The continuity equation and the energy-momentum conservation law have the form 

" ^' Ji^) = '^^^'^ +Y^kw'^k)^ Jf^-^eW", (1) 

r.7 =0, T^'^'^{P + e) u^u'' + Pg^"" + Y,k (wf,) w^^^ + w^^^ u"" + /ifc w^^ u^) . (2) 

Here and below the subscripts i and fc refer to nucleons: i, k = n, p; the subscript / runs over all particle species, 
I = p, e. Unless otherwise stated the summation is assumed over the repeated nucleon indices i, k and over the 
spacetime indices fi, and a. In Eqs. ([T]) and (HJ, j^^j^ is the particle four-current for species Z; T^'' is the energy- 
momentum tensor; g^'^ is the metric tensor; ni and /i; are the number density and chemical potential of particle 
species I, respectively; e and P are the energy density and pressure, respectively. The four-vector is related to 
the superfluid four- velocity u^^^^ by equality: w^^-^ = fii(v'^^^^ — u'*). In what follows, all the thermodynamic quantities 
(e.g., ni and s) are deflned (measured) in the reference frame comoving with the normal component of the fluid [in 
this frame = (1,0,0,0); notice that is normalized in such a way that m^m'' = —1]. This imposes an additional 
constraint on the four- vector w'^-s^ [III illl i 

Wpwf,) = 0. (3) 

Using this constraint and Eqs. ((IJ and ^ one immediately finds that, in a relativistically invariant form, the particle 
number densities are defined as ni = ^u^j'^^j 35], while the energy density is e = u^u^T^^'^ . 

The matrix Yik in Eqs . ([IJ and ([5]), named the relativistic entrainment matrix, is a generalization of the concept 
of superfluid density [l5| to the case of relativistic mixtures. In the nonrelativistic theory, a similar matrix was first 
considered by Andreev and Bashkin [s^. The matrix is expressed in terms of the Landau parameters Fl^ of 
asymmetric nuclear matter and universal functions of temperature, <l>i, as described in [32j . Generally, it can be 
presented as a function of density p and the combinations T /Ten and T /Tcp'. Yik = Yik{p, T /Ten, T /Tcp), where Tcp is 
the proton critical temperature. 

For illustration, in Fig. 1 we plot the elements Yik as functions of temperature T for the fixed baryon number 
density rn, = 3no = 0.48 fm~'^. The matrix elements are normalized to the constant Y — 3no//in(3no) ~ 2.60 x 
10**^ erg~^ cm~'^, where no = 0.16 fm~'^ is the density in atomic nuclei and /i„(3no) is the neutron chemical potential. 
Here and below, in all our calculations we use the parameterization [s^j of equation of state (EOS) of Akmal et al. [11] 
(hereafter APR). The vertical dot-dashed lines indicate the critical temperatures for neutrons, T^n = 6 x 10® K, and 
protons, Tcp = 5 x 10^ K. Let us stress that the actual values of baryon critical temperatures as well as their density 
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FIG. 1: Normalized symmetric matrix Yi^jY as a function of T for rth = 3no = 0.48 fm~'^. The normalization constant 
Y = 3no/Atn(3no) ~ 2.60 x 10^^ erg~^ cm"''. Solid lines show the elements Yi^jY {i,k = n,p). Vertical dot-dashed lines 
indicate the baryon critical temperatures Ten = 6 x 10* K and Tcp = 5 x 10^ K. 



profiles are poorly known. Different authors obtain different results for T^n and Tcp depending on what microphysical 
input they use and how they take into account the many body effects [3, Q . Because of a large uncertainty in the 
models of baryon critical temperatures, in this paper we will treat them as free parameters. Our models of Ten and 
Tcp do not contradict these, available in the literature. 

One sees from Fig. 1 that at low temperatures (T < 10^ K), the elements of the matrix Yik are nearly independent 
of temperature and coincide with their values at T = 0. With increasing temperature, they start to decrease and 
vanish after the corresponding critical temperature is reached. Thus, the matrix Yik strongly depends on T near the 
superfluid phase transition. 

Along with the continuity equations and the energy-momentum conservation, one should impose an additional 
constraint on the superfluid velocities, following from the potentiality of superfluid motion 



9" 



-I- q.A'' + muf" 



(4) 



where qi is the charge of a species i, and is the four-potential of the electromagnetic field. 

To close the system, the hydrodynamic equations should be supplemented with the second law of thermodynamics, 
which has the form [ll|, 



Yik 

de = T dS + He dric + Hi drii H — — d 



(5) 



where S is the entropy density. 



B. Radial pulsations 

The metric of a radially pulsating nonrotating symmetric NS has the form (see, e.g., [39j): 

ds2 = -e'^dt^ + r^dn^ + e^ dr^ , (6) 

where r and t are the radial and time coordinates, respectively; dil is the element of a solid angle in a spherical 
coordinate system with the origin at the stellar center. The metric coefficients v{r,t) and A(r, i) depend only on r 
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and t. In what follows, the quantities related to a star in the hydrostatic equilibrium will be denoted by the subscript 
'0'. For example, the metric coefficients of an unperturbed star are designated as voir) and Ao(r). 

Using Eqs. (l)-(5) the authors of GA06 derived the system of linear equations that describes radial pulsations of a 
general relativistic superfluid NS (all perturbations were assumed to be small and dependent on time t as exp(ia;i), 
where ui is the pulsation frequency): 

e^"-''" uj^ mofino^b = ^ + ^ (^^0 + + 
2'^' d7"2'"'°^"°^^ I d7+rj d;:(^o + ^o)' 

(8) 



Here, the quantities fc, dP, dfin, and (5i^ are, respectively, perturbations of e, P, fin and the metric coefficient v{r,t)] 
^ is the Lagrangian displacement of the normal liquid component, related to its velocity by [s^ 

ui=e-''°/2||=ic^e-'^«/2e (9) 

Furthermore, Zi is the Lagrangian displacement of a superfluid component i relative to the normal component, defined 
as 

and the function ^j, is proportional to the baryon current density j^J^^ = nt,u^ + Ynkw'^y^ + [^^'^ Eq. ([1])], and 

is given by 

tb = 4H \ . (11) 

nbo Ubo 

Equation ([7|) is the relativistic Euler equation generalized to the case of superfluid stars. A similar equation has 
been derived for normal matter by Chandrasekhar [s^. Equation (|5]), specific to superfiuid matter, is a consequence 
of the potentiality of superfluid motion [see condition (|4])]. 

Because we are interested in oscillations whose frequencies are substantially smaller than the electron plasma 
frequency, one can assume that the quasineutrality condition Hp — Uf. is always satisfied. As a consequence, the 
electron and proton currents must be equal, j'^^-^ = j'^'^-j. Using now expression ^ and the definition (llOp . one obtains 
an additional constraint that relates z„ and Zp, 

A^fco Ypk Zk = 0. (12) 

In this study, it will be convenient to transform Eq. ([8|) and rewrite it in a somewhat different form. For this 
purpose, we multiply Eq. ^ by nf,o exp(— z/o/2) and subtract Eq. JT]) from the result. After some algebra one gets, 
using equation (54) for dSv/dr from GA06 (see also [H, [slj), 

\o-ua 2 ( f^koYnkZk\ dSfi 

Aino n-bo e " ° u [ z„ = rieo ^— , (13) 

V "60 / or 

where = ^„ — /Ltp — fig is the disbalance of chemical potentials. Solving the system of equations Q, ([H]), and (|13p . 
we can find the eigenfunctions ^, z„, and Zp, as well as the spectrum of eigenfrequencies uj. In order to do so, it is, 
however, necessary to express the quantities Se, SP, and (5/i in terms of the functions Zn, and Zp and to formulate 
the boundary conditions for Eqs. ([7]) and (|13p . 

Owing to the quasineutrality condition, = Hp, any thermodynamic quantity can be presented as a function of 
two variables, say, nb and Tig (the dependence on the quadratically small scalars w'^{-fW(k)fi and on T can be neglected). 
Because we consider small-amplitude pulsations, the pressure perturbation and the chemical potential disbalance can 
be expanded in Taylor series near their equilibrium values, 

dP{nbQ,neo) ^ , dP{nbQ,neo) . 
dP = dub-i One, (14) 

dUbo OUeO 

O/i = OHb-i One- (15) 

OUbO oneo 
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Similarly, in view of the beta-equilibrium condition S^iq = /i„o ~ MpO ^ Meo = 0, valid for an unperturbed star, one has 
from Eq. (O (see also GA06), 

5e = finoSnb. (16) 

Using continuity equations ([T]), the variations Sni, and Jjig appearing in Eqs. (|14l) - (jl6p can be expressed in terms of 
the Lagrangian displacements ^, z„, and Zp, 



- ~^d-,{r'-^o^,e-^"/'), (17) 

j,2 Qj, 



fee = -^|:(r2neofe-''°/2). (18) 



C. Boundary conditions 



Prior to formulating the boundary conditions for Eqs. ([7]) and ((T5)) . we should refine the model of a star whose 
pulsations we analyze. We consider a two-layer model constituted by two regions: an inner region, in which neutrons 
are superfluid, and an outer region, in which neutrons are normal. In the inner region (for which r < Rq), the 
functions Zi are nonzero; we name this region 'superfluid'. In the outer region (for which r > Rq), all matter 
components move with the same velocity u^^, i.e., the fluid demonstrates a nonsuperfluid behavior [even if protons 
remain superfluid, it follows from (|T2t that Zp — 0]; we name such a region 'normal'. In addition, we assume that 
the thermodynamic functions (e.g., the particle number densities ni) do not experience any abrupt change at the 
superfluid-normal interface. Then one can formulate the following boundary conditions for such a stellar model (see 
also GA06). 

(1) Vanishing of the pressure at the stellar surface. 

= 0, (19) 

r=R 

where R is the circumferential radius of the star. 

(2) Existence of all the derivatives with respect to r at the stellar center, which means that the following limits 
must be finite 



dr 



lim ^t/r < oo, lim fiioYni Zi/r < oo. (20) 

r— )-0 r— >-0 

(3) Continuity of the energy current density across the superfluid-normal interface. As shown in GA06, this current 
is directed along the radial coordinate and is proportional to rif,oA*rioC + A'jioMio^ni-Zi = fiboA*?i0^b- Hence, because of 
the continuity of the thermodynamic functions, one has 

(4) Continuity of the electron current density across the superfluid-normal interface, ^(-Ro ^ 0) = C(^^o + 0), which, 
combined with Eqs. (HH), ([H]), and (HI]), leads to 

^J■ioYni Zi |r=_Ro= 0. (22) 

The solution to the system of equations ([7]), (IT2]) . and (|13p must satisfy these four boundary conditions. 



III. PULSATION SPECTRA 



Using Eqs. ([7]), (IT2t . and ([13]), we determined the spectrum of eigenfrequencies w of a radially pulsating NS as a 
function of its internal temperature T. In our numerical calculations, we considered NS with the mass M = 1.4M0, 
circumferential radius R = 12.17 km, and central density pc = 9.26 x 10^'' g cm~'^. We used APR EOS in the core, 
and Negele and Vautherin EOS in the crust [10]. For the chosen model of a star, the core-crust interface lies at a 
distance i?cc = 10.88 km from the stellar center. Following GA06, free neutrons in the inner crust were assumed to 
be normal, while in the core both nucleon species were allowed to be superfluid. 

To simplify calculations, the red-shifted critical temperature of protons was chosen to be constant throughout the 
core (i.e., independent of the density p) and equal to = Tcpe'^°^^ = 5 x 10^ K. Note that, the pulsation spectrum 
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only weakly depends on a particular form of the function T^{p) and even on whether or not protons are superfluid. 
The point is the only quantity in superfluid hydrodynamics, which depends on the critical temperatures T^n and T^p, is 
the entrainment matrix Yik ■ The element y„„ of the matrix is mostly determined by the neutron critical temperature 
Ten and is nearly insensitive to Tcp [32j. By contrast, the elements Ypp and Ynp strongly vary with T^p, but the terms, 
depending on these matrix elements, are small in pulsation equations and can be neglected to a good approximation. 

To illustrate this statement let us notice that, as follows from Eqs. ([7]) and (fT3|) . Ynp appears there only in the 
combination /i„ y„„ z„ + /Zp Ynp Zp (and Ypp does not explicitly enter these equations) . From the quasineutrality 
condition (|12p . we obtain Zp = —fJ.nYnpZn/{fJ-pYpp), and, therefore, the second term in this combination differs from 
the first one by a factor Y^p / (YnnYpp) . As follows from Fig. 1, Yj^p/{YnnYpp) ^ 0.01 and, consequently, this term can 
be indeed neglected. 

The red-shifted temperature of unperturbed star, T°° — Te'^°^^, was assumed to be finite and constant throughout 
the stellar core. The latter is a necessary condition for a superfluid star to be in the hydrostatic equilibrium [see 
equation (30) of GA06]. The influence of finite internal temperature on the pulsations of NSs is twofold. First, it 
determines the entrainment matrix Yik{T) (see Fig. 1). Second, it regulates the size Ro of the superfluid region, for 
which T < Tcn{p)- the higher the temperature, the smaller the region (see Fig. 2). 

We considered four models showing different behavior of Ten (p) and allowing for one or both of these temperature 
effects. 

Model 1. The red-shifted neutron critical temperature T^ is density independent and equals T^^ — Tcne'^°^^ = 
6 X 10* K. This model assumes that superfluid region coincides with the stellar core at T°° < T^. Within this region 
we use the flnite- temperature hydrodynamics, i.e., take into account the full dependence of Yik on T. This model is 
identical to the model considered in GA06. The only difference is that the authors of GA06 used an approximate 
expression for the matrix Yik [see their equation (17)], whereas in the present study we use results of a self-consistent 
relativistic calculation, performed in 's^. However, one can check that the expression for Yik from GA06 is a good 
approximation for calculation of the pulsation spectrum. 

Model 2. T^ depends on density and has a maximum at the stellar center, T^{r = 0) = 6 x 10* K. The size of 
the superfluid region is temperature dependent and the boundary Rq of this region is determined by the condition 

r,^(i?o) = r-. (23) 

However, within the superfluid region we use the zero-temperature hydrodynamics, that is, the entrainment matrix 
in this model does not depend on T, Yik = YikiT = 0). This model has been used in a number of studies f23l. [28l. [29j . 



In particular, in Ref. [23[ it was claimed that the model adequately accounts for the main temperature effects. 

Model 3. The same as model 2, but within the superfluid region we use a hydrodynamics which is valid at flnite 
temperatures; Yik is temperature dependent. 

Model 4. The same as model 3, but with having a maximum at a density exceeding the central density pc- At 
the stellar center, T;^{0) = 6 x 10* K. 

It should be stressed that, for each of these models, the results obtained below in this paper will remain qualitatively 
the same if we choose another neutron critical temperature T^{0) in the stellar center or if we vary the slope of the 
density profile T^(p) (for models 2, 3, and 4). The only important thing that discriminates between the models 
and that should not be altered is the condition [dT^{p) /dp]\r=o = for models 2 and 3 and [dT^{p)/dp]\r=o > 
for model 4. As it will be demonstrated in Sees. HI and IV, this condition infiuences the pulsation spectra and the 
asymptotes for eigenfrequencies at T°° T^(0). 

For illustration. Fig. 2 shows the red-shifted T^ (dashed curve) and the local Ten (solid curve) neutron critical 
temperatures as function of density p (left panel for models 2 and 3; right panel for model 4). Dot-dashed lines 
indicate the central density of the star and the density at the crust-core interface. Notice that, for all the models, a 
sufficiently cold star (at T°° < 6 x 10* K) consists of superfluid inner and normal outer regions. 

In Fig. 3, for each of the four models, we present an eigenfrequency spectrum of a radially pulsating superfluid NS. 
The alternate solid and dashed lines show eigenfrequencies u) (in units of ujq = c/R) as functions of T°° for the flrst 
six pulsation modes. The dot-dashed line indicates the red-shifted neutron critical temperature at the stellar center, 
T^{0) = 6 X 10* K. The dotted lines turning into thin solid lines at a temperature T°° > T^{0) show first three 
radial oscillation modes of a normal star (I, II, and HI). The spectrum was not plotted in the shaded region. 

It can be seen from the figure that at < T^{0) there exist two types of oscillations. Oscillation frequencies 
of the first type are nearly temperature independent and well coincide with the corresponding eigenfrequencies of a 
nonsuperfluid star of the same mass (we name this type of modes 'normal'). Oscillation frequencies of the second 
type strongly vary with temperature (for T > 5 x 10^ — 10* K); the corresponding modes are named 'superfluid'. 
With growing T, neighboring modes become closer and closer until they form avoided crossing. Near such avoided 
crossings any superfluid mode turns into a normal mode and vice versa. At low enough temperatures (T < 5 x 10^ 
K for models 2, 3, and 4, and T < 10* K for model 1), the eigenfrequencies cease to depend on T and reach their 
asymptotic values. The reason for that is, flrst, the matrix Yik becomes temperature independent at T°° < 10* K 
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and second, the radius of the superfluid region reaches its maximum value (corresponding to the core-crust interface) 
at T°° = 5 X lO'' K and stops to vary as T is lowered further. 

For the superfluid type of pulsations, the decrease in the size Rq of the superfluid region with increasing T°° leads 
to a rise of the eigenfrequencies uj (this is clearly seen for model 2 and stems from the fact that, as shown below, 
uj cx I/Rq). At the same time, Ynk — > with increasing T (Fig. 1), which results in decreasing of uj (this is exemplified 
by model 1 and follows from the fact that uj oc V^rm at T°° — > T^^(O), see below). In the models 3 and 4, these 
two effects compete. Asymptotes for the eigenfrequencies at T°° — T^{0) will be studied in detail in Sec. IVB. At 
temperatures higher than T^{0), the star pulsates as a normal one [even if protons still remain superfluid, all the 
liquid components move together by virtue of the quasineutrality condition (|12p ]. 

As follows from Fig. 3, the pulsation spectra differ substantially for different models. For example, for most realistic 
models 3 and 4 in the limit T°° — > r^^(O), the eigenfrequencies either vanish (model 3), or decrease, approaching some 
certain finite values (model 4). The eigenfrequencies for model 1 behave similar to those for model 3. At the same 
time, the eigenfrequencies for model 2 behave in a drastically different way. At T°° — ?> ^^^(O) they grow rather than 
decrease and tend to infinity as Rq 0. This example clearly demonstrates that it is insufficient just to vary the size 
of the superfluid region with T, using the zero-temperature hydrodynamics (as in model 2). Such an approach gives 
qualitatively incorrect spectra. Thus, it is of principal importance to take into account the temperature dependence 
of the entrainment matrix. 



IV. EIGENFREQUENCIES AT r°° T^(0) 
A. Approximate splitting of pulsation equations 

As shown in Refs. ^33*, '33], the equations that govern global pulsations of superfluid NSs can generally be represented 
as two weakly coupled systems of equations. The coupling parameter s for these systems is given by 

^ _ ?ieO dP{nbo,neo)/dneo ^^^^ 
ribQ dP{nbo,neQ)/dnbo 

and is small for realistic equations of state, \s\ ~ 0.01 — 0.05. In particular, s for APR EOS is shown in Fig. 4. In 
the case of radial pulsations, the smallness of the parameter s means that the pulsation equations ([7]) and can 
be treated as independent to a very good approximation (see [H, [13] )• If we assume s = 0, then Eq. ([7]) can be 
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FIG. 3: Eigenfrequencies a; (in units of ixjq = c/i?) versus Tg" — r°°/10*, calculated for models 1, 2, 3, and 4 (see the text). 
Solid and dashed curves show the first six pulsation modes. The red-shifted neutron critical temperature Tc^(O) in the stellar 
center is indicated by the vertical dot-dashed lines; the horizontal dotted lines turning into solid lines at T°° > rc^(O) show 
the first three eigenfrequencies (I, II, III) of a normal star of the same mass (notice that, for model 2 we show the first six 
eigenfrequencies of a normal star). No spectrum was plotted in the shaded region. 



separated from Eq. p3p and written only in terms of the variable ^\,: 

e " " uj nfcoAi«o6 = h oPnorm 3- ;rAo + fo 

or dr \2 



where 



1 di>o 1 / di'o 1 \ d 

^ 2 dT ~ 2 "^"^"0^^ i d7 + r j + ''o) ' 



xrj _ dP{nbo,neo) aP(n6o,n-eo) . 

O-Tnorm = T{ (""b H (""e 

OUbO OUeO 



e 



j,2 Qj, 



|:(.^n.o6e--/^). (26) 
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FIG. 4: Coupling parameter s versus p for APR EOS. The vertical dotted lines indicate the central stellar density and the 
crust-core interface. 



In this form, the equation fully coincides with that describing oscillations of a nonsuperfluid star. Its solutions are 
eigenfrequencies w and eigenfunctions of a normal star. Thus, Eq. ((25|) approximately describes the normal modes 
of oscillations in our problem. As for the superfluid modes, they are described by an equation depending only on 
'superfluid' Lagrangian displacements z„ and Zp (ssl. l33|. 

Mno"boe " "UJ Zn = neo — , (27) 

V nbo J or 



where 



dfisa = -t: 5— -5- fiioYniZiC °' 

dUe dr \ ribo ) 



(28) 



Eq. (P71) can be derived from ([T^ if we put ^6 = in the latter. Indeed, if some frequency is an eigenfrequency for 
the system of equations (|13p and (|25p . but is not an eigenfrequency for Eq. ([25]) . then this is only possible if f;,, the 
eigenfunction of Eq. (|25p . identically vanishes. (This consideration is strictly valid when s = 0. For small but finite s 
one has < ^ w -[liaY-ai Zi/ribQ.) 

The approximation of independent Eqs. (1251) and (|27p well reproduces (with an accuracy of 1.5 — 2%) the frequency 
spectrum of a pulsating NS. However, this approximation does not describe avoided crossings of neighboring modes 
(see Fig. 3). The reason is near these points Eqs. ([7]) and (IT3|) are strongly interacting [S^Hj], so that approximate 
treatment based on the decoupled Eqs. (1^ and (P7|) is unjustified; instead of avoided crossings, Eqs. (pSj) and (P7)) 
predict crossings of modes. It should be noted that the existence of two types of modes, superfluid and normal. 



has been noted previously by a number of authors (see, e.g., Refs. [IJ, l28|, |41[ ) . However, an adequate analytical 
interpretation this fact first received in Ref. [33| . 



B. Asymptotes for eigenfrequencies 

Using the s = approximation described above, let us consider the asymptotic behavior of eigenfrequencies at 
r°° — ;> T^[()). For all of the four models, any given mode experiences superfluid- type pulsations near the neutron 
critical temperature and, therefore, to derive the asymptotes it is necessary to use Eq. (P7)) . 

As the temperature approaches the critical value Tc„, the elements Yni{p,T /Tcn,T /Tcp) of the entrainment matrix 
tend to vanish according to the law [S^Ull: Yni — Fni (1 — T/Tcn), where Fni is some (known) function of p, T/Tcn, 
and T/Tcp, which are in turn the functions of radial coordinate r. Taking into account that T /Td = T°° /T^ , the 
asymptote for Yni can be rewritten as 

Yn^ = F„^{r,T°°){l - T°°/T^) 0. (29) 
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Owing to the smallness of the matrix elements Y„i near the neutron critical temperature, one can neglect the second 
term in the left-hand side of Eq. (|27p when analyzing asymptotes for models 1, 3, and 4. 

In addition, the size Rq of the superfiuid region in models 2, 3, and 4 decreases with increasing temperature and 
vanishes at T°° = T^{0). The limiting behavior Ro{T°°) can be found from the condition (1^^ and is determined 
by the profile of the neutron critical temperature. Expanding T^{p) in the Taylor series in the vicinity of the stellar 
center and taking into account that for models 2 and 3 [dT^/dp]\r=o = while for model 4 [dT^^/dp]|r=o 7^ 0, one 
gets, for models 2 and 3, 



(30) 



and for model 4, 



?^™M-r,^(0)«p(r)-p,(xr2 



(31) 



To obtain the asymptotes (PH)) and (PT|) we used the fact that in the stellar center [dp/dr]|r=o = 0. Now, using Eqs. 
([30)1 and dST]) and the condition one finds how Rq scales with T°° at T°° T^^iO). For models 2 and 3 



i?o(x[l-r-/T-(0)] 



1/4 



(32) 



while for model 4 



i?ocx[l-^-/T-(0)]l/^ 



(33) 



Making use of Eqs. (I12)-(ISS1), it is possible to analyze the dependence w(T°°) at T°° T^{0) for different models. 
We start with the model 1. Then, the only quantities in Eq. (|27p that depend on T are the matrix elements Yni and 
the frequency w (we remind that the size Rq of the superfiuid region is fixed for model 1). Using Eq. and the 
fact that T°°/T^ does not depend on r for this model, one can factor the temperature dependence of Yni outside the 
differentiation sign on the right-hand side of Eq. (|27p . The result is 



one V or \ Ubo 



(34) 



The function Fm{r,T°°) « Fm{r,T^) in the immediate vicinity of 
obtains from Eq. (p4|) . 



■oo 

C7l ■ 



Taking this into account, one immediately 



w oc v^l - T^/T^7 ^ 0, T° 



(35) 



Now let us inspect model 2. This model assumes that the matrix Yi^ is temperature independent, while i?o 
decreases with increasing T . In that case, the asymptote for the function uj(T°°) can be determined by introducing a 
new dimensionless variable, r = t/Rq, which varies within the range from to 1. Rewriting Eq. (|27p employing the 
new variable, one finds [notice that we retain the second term on the left-hand side of Eq. (P7)) for this model] 



Pnoribo e 



Rq 



PkoYnk Zk 



"60 



d_ 

df 



zir- — PioYn 

One or \ nto 



(36) 



At T°° sufficiently close to T^^(O), when the size Rq of the superfiuid region is rather small, such functions as //iOj 
life, nio, and vq are almost independent of the radius r [so that one can put, e.g., p.io{r) /iio(O)]. Thus, the only 
quantities varying with temperature in Eq. (|36p are uj and Rq, and, therefore, the asymptote takes the form 



1 

UJ cx — oc 

Rq 



1 



[1 - r-/T-(o)] 



1/4 



OO, T-^r-(o). 



(37) 



To obtain this asymptote we also used Eq. ([32|). 

Finally, let us consider models 3 and 4. In these models, both Yik and Rq depend on T°°. So, we need to combine 
our reasoning for the first and second models. The pulsation equation for models 3 and 4 is given by the same 
Eq. p6l) but without the second term oc Yni on its left-hand side; as it is argued above in this section, this term 
is small at T°° — and can be omitted. Again, as for model 2, in Eq. (|36p one can replace such quantities 
as iiio{r) by their values /iio(O) in the stellar center. The same, of course, cannot be done for the matrix element 
r„i(r,T°°) ^ Fm{r,T°°) [1 - T°° /T^{r)\, because it changes substantially in the superfiuid region. Using Eqs. ^ 
and l|29p - ((5T|) it is easy to verify that y„i(r, r°°) can generally be presented as 



Y 



F™(0,T-(0)) [1 - T-/T~(0)] (1 - f*) + 0([1 - T^/T^mn 



(38) 
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for model 3 and as 

= i^„(o,T-(o)) [1 - r-/T-(o)] (1 - + o([i - r-/r-(o)]2) (39) 

for model 4. The temperature dependence of can now be easily factored out of the differentiation sign on the 
right-hand side of Eq. (j36|) . In this way, one obtains the following asymptote for w, 



i?0 



'^^(i-T?;^) ->0' r-->r-(o), (4i) 



Now, employing Eqs. (|32|) and (|33|) , one finds for model 3 
and for model 4 

finite, r°°^T^^(0). (42) 

It should be emphasized that, owing to a weak dependence of the function F„i(r, T°°) on T°°, the eigenfrequencies 
for the model 4 do not tend to a constant, but remain weakly dependent on temperature. One can verify that the 
asymptotes obtained for all the models well describe our numerical results presented in Fig. 3. 

In this study, we only consider profiles {p) that have a maximum at densities higher than, or equal to the central 
stellar density pc- This enables us to consider only two regions in the star: inner superfluid and outer normal. However, 
a situation in which the maximum lies at p < pc can easily be imagined. In this case, there will be generally three 
regions in a star: inner and outer normal regions and a superfluid region in between. The size i?o of the superfluid 

1/2 

region will decrease with increasing temperature in the same manner as in model 4, i?o oc (1 — /T^^^^) ' (where 
^mmax is the value of at the maximum). Therefore, at T°° — > T'mmax the eigenfrequencies for such star will have 
the same asymptotes as for the model 4, w — > finite. The spectrum for a three-layer star is analyzed in Ref. [43j . 



V. EIGENFUNCTIONS 



In this section we compare the main properties of superfluid and normal pulsation modes. For that we analyze the 
eigenfunctions for both types of modes. For brevity, we consider only the third model. 

Figure 5 shows the eigenfunctions ^, z„, Zp, and ^j, (normalized to the stellar radius i?), calculated at T°° = 
1.2 X 10^ K versus r (in units of i?). The oscillation amplitude was chosen in such a way that the energy of radial 
pulsations is -Epuis = 10^^ erg [see Eq. (|44l) for the definition of iSpuis]. For a fundamental mode of a normal star this 
energy corresponds to an oscillation amplitude ^ ^ lO^'^r. It could be also interesting to note that this energy is of 
the order of the thermal energy of a normal star with = 10^ K [S]. Each mode (I,...,VI) in Fig. 5 is plotted in a 
separate panel. The boundary i?o between the superfiuid and normal regions, and the crust-core boundary i?cc) are 
indicated by vertical dotted lines. The solid line corresponds to the Lagrangian displacement ^ of the nonsuperfluid 
component; the dashed and dot-dashed lines correspond to relative Lagrangian displacements z„ and Zp, respectively; 
the heavy solid line - to the Lagrangian displacement ^h, which is proportional to the baryon current density. We do 
not plot Zn in the nonsuperfluid region of the star since all neutrons are unpaired there. We also do not plot Zp in the 
crust, because it contains no free protons. The quantity in the normal region of the star coincides with ^. 

At the chosen temperature (r°° = 1.2 x 10* K), the modes I, III, IV, and VI experience oscillations of the superfluid 
type and are approximately described by Eq. (j27p . The modes II and V are of normal type and are approximately 
described by Eq. (j25[) : their eignefrequencies roughly coincide with the flrst and second modes of a normal star of 
the same mass. It can be seen that, as expected, the baryon current in superfluid oscillations, characterized by the 
quantity ^f,, is almost zero. 

It is noteworthy that the crust of a star, experiencing superfluid oscillations, remains nearly at rest, i.e., pulsations 
are mostly localized in the superfluid region. This point is illustrated in Fig. 6, which presents the ratio of average 
oscillation energy densities in the nonsuperfluid enorm and superfluid Egfl regions as a function of T°° for the flrst six 
pulsation modes (I,..., VI). Modes I, II, and III are shown, respectively, by the solid, dotted, and dashed lines in the 
left panel; the other modes are shown in the right panel. The thin solid curve is plotted for the first two modes of a 
normal star (Infh and Ilnfh). The latter modes depend on temperature only because of the variation of i?o with r°°, 
which is artificially assumed to be the same as for a superfluid star. 

It can be seen from Fig. 6 that, for normal-type oscillations, the ratio enorm/ssfl coincides very well with the result 
for normal star. However, immediately after the star starts to experience pulsations of the superfluid type, Enorm/^sfl 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 

1^ y ly^ 



FIG. 5: Eigenfunctions ^, 2„, Zp, and ^i, (normalized so that Spuis ~ 10** erg) of a superfluid NS at T°° = 1.2 • 10* K versus 
r (in units of R), calculated for the first six pulsation modes (I,..., VI) of the model 3. Solid curves correspond to ^, dashed - 
to Zn, dot-dashed - to Zp, heavy solid lines show ^6. Vertical dots indicate the crust-core boundary and the superfluid-normal 
interface. 

falls nearly to zero, so that almost all pulsation energy will be localized within the superfluid region. Thus, we see 
that the superfluid pulsation modes practically do not appear on the stellar surface and therefore it should be difficult 
to observe them by detecting the modulation of electromagnetic radiation from the star. 



VI. DAMPING TIMES 



In this section we calculate characteristic damping times for radial pulsations due to the shear viscosity and nonequi- 
librium reactions of mutual particle transformations. It is well known that such reactions generate an effective bulk 
viscosity (see, e.g., [Ulli^ and references therein). 

The e-folding time is defined as 



{E-pu\s) 



(43) 



where (iJpuis) is the damping rate for Spuis, averaged over the period 27r/aj. Here and below, the angular brackets 
denote averaging over the pulsation period. 

The expression for the pulsation energy of a superfluid symmetric nonrotating star has the form 



1 



MiOA^fcO Yik Zi(a)Zkia) + MiOMfeO ^ik S,{a) (Zi (a) + ^fe(a)) + (^0 + ^o) S,fa) 



(44) 



where 2:^(0), ^fe(o)j and ^(a) are the amplitudes of Zi, Zk, and ^, respectively. Using the symmetry of the matrix Yik, 
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FIG. 6: The ratio of the pulsation energy densities in the nonsuperfluid and superfluid regions of a star versus Tg° , calculated 
for model 3 (see the text). The left panel shows the first three modes (I, II, and III), the right panel shows the modes IV, V, 
and VI. Solid thin lines present the same ratio for the first two pulsation modes of a normal star (Infh and Ilnth). The redshifted 
neutron critical temperature in the center of the star is indicated by the vertical dot-dashed lines. 



the quasineutrality condition (|12p . and the definition (jlip of ^j,, we can simplify the expression for -Epuis, 

EpvtXs ^ 2 J |.M«OMfcO Ynk Zk (a) [ Zn (a) — j +nbO^J■nO 4 (q) 

x47rr2e3^«/2-^o/2^^^ (45) 

For a stellar model considered here, the direct Urea process is forbidden, and the only effective processes of particle 
transformations are the neutron and proton branches of the modified Urea process, 

+ >p + n + e^+I7e, p + n + ^ n + n + i/e', (46) 
n + p ^ p + p + +Ve, p + p + ^ n + p + Ue- (47) 

These reactions lead to dissipation of iJpuis and to heating of the star. Assuming that pulsations are subthermal ((5/i <C 
fcT), the conversion rate of the pulsation energy into the thermal energy (per unit volume) via the nonequilibrium 
reactions (j46|) and (gll), is given bv [l7l[3ll43| 

Qbuik = AJ^^^ (48) 

where the chemical potential imbalance 5/i can be expressed in terms of the above-calculated Lagrangian displacements 
4 and Zp with the help of Eqs. ([T5|) . p7)) . and and A = Ai -I- A2 is the sum of the reaction rates and (|47t . 
The reaction rates and the corresponding reduction factors were calculated for the modified Urea process in Ref. [4^] 
(see also Ref. [l^l , in which a correction was introduced to the reaction rate for the proton branch) . 

The energy converted into heat per unit volume per second due to the shear viscosity is given by the same formula 
as for normal matter [43 |. 

Q..... = \^^'-^ {~^4 + 2^ + r^^-p-)\ (49) 
6 r'^ \ or ar J 

where 77 is the shear viscosity coefficient. In our numerical calculations we used 77 from [i^. Because the shear 
viscosity is mostly determined by the electron-electron and electron-proton collisions, we neglect the contribution of 
neutron collisions and use equation (37) from ^48] when calculating 77. In addition, we use equation (83) from the 
same reference to take into account the reduction of the particle collision frequency (and, consequently, the increase 
in the shear viscosity) by the proton superfiuidity. 
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FIG. 7: Damping time r versus T°° for the first six pulsation modes of model 3 (see the text). The thin solid lines indicate 
r for the first two modes of a normal star (Inm and Ilnfh). Because of numerical problems t was not calculated in the filled 
region. Other notations are the same as in Fig. 6. 



To obtain the dissipation rate of the pulsation energy, one has to average Eqs. (|48| and (|49| over the pulsation 
period and then to integrate them over the stellar volume [44l] . 

(-Epuls) =- j ((Qshcar) + {Qx,nl\.)V'' AtTT^ dr, (50) 

Having thus expressed all the quantities in Eq. (|43p in terms of the already calculated eigenfunctions, we can present 
numerical results for the e-folding times. As in Sec. V, we consider only the model 3. 

Figure 7 presents t as a function of T°° , for the first six pulsation modes (three modes in each panel) . As in Fig. 
6, the modes are shown by heavy soUd, dotted, and dashed lines. Thin solid lines show r calculated for the first two 
modes Infh and Ilnfh by using normal-fluid (ordinary) hydrodynamics. In these calculations, superfluidity was taken 
into account only at calculating the kinetic coefficients (such an approximation is frequently employed in the literature 
to estimate the effect of superfluidity on the dissipative processes in NSs). The damping times were not calculated for 
the modes IV, V, and VI in the immediate vicinity of T^^(O) because of problems with numerical calculations (this 
region is filled grey). 

It follows from Figs. 3 and 7 that t for normal-type oscillations may differ severalfold from that calculated using 
the normal-fluid hydrodynamics. The difference is stronger at lower temperatures. This fact may seem surprising 
if we recall that the eigenfrequencies and eigenfunctions ^^(r) approximately coincide for normal modes and for 
corresponding modes of a nonsuperfluid star. To explain the fact, one should notice that the function ^(r) generally 
differs in superfluid and normal stars. In additon, the relative motion of superfluid components with respect to the 
normal one {zi ^ 0) also modifies r. 

Superfluid pulsation modes are damped faster than the normal modes, so that r drops by 1 — 3 orders of magnitude 
each time the star switches (e.g., in the course of thermal evolution) from normal- to superfluid- type regime of 
pulsations (see Figs. 3 and 7). To explain such a fast damping, let us note that it occurs (for r°° < 6 x 10^ K) mainly 
due to the electron shear viscosity. For a nonsuperfluid matter the electron shear viscosity coefficient r]^ is smaller in 
the outer stellar layers than in the inner layers. Proton superfluidity further increases ?7e in the inner layers of the 
star. Now, if one looks at Fig. [51 one will notice that superfluid modes are localized in the superfluid core, while the 
amplitude of normal modes (modes II and V in that figure) grows away from the stellar center (with the maximum 
in the crust, where "qe is small). In view of these facts the reason for the difference of damping times becomes clear; 
quantitatively, this result follows from a careful examination of Eqs. P5t . and ([50]) . 

In order to understand the limiting behavior of r at high temperatures [T°° — >■ T|,^(0)], let us analyze its asymp- 
totes for different models. As already discussed in Sec. IV, any given mode near the neutron critical temperature 
experiences superfluid- type pulsations, approximately described by Eq. (j27p . For this type of pulsations ^f, is small. 
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<^ Mfeo Ynk Zk/nb, SO that one derives from Eq. (|45t the following estimate for -Epuis at T°° — )• T^{0) 

MzO -^^i -^■i (a) 



I (a) 



?^fcO 



Rl (51) 



Here and below, we neglect the factors of the order of unity in aU formulas and replace integration over the stellar 
volume with the factor Rq] the quantities z„(a) and 2p(a) in Eq. (jSTj) denote Zn(a) and Zp(Q), averaged over the 
superfluid region. Taking into account that for superfluid modes ^ ~ —fJ'koYnkZk/nbo one obtains, using Eqs. (pS)) 
and (j49p . the following estimate for the dissipation rate due to, respectively, the shear viscosity and nonequilibrium 
modified Urea reactions (bulk viscosity) 



/ MfcO Ynk Zk (a) 



2 



(gpuis sh) - -w' ( — -"^"^ ) i?o, (52) 



An^„f^^^^i^^^Vf|^Vi?o. (53) 



{Epu 



Is b; 



In view of Eqs. ((5T1) - ([55|) one finds for the corresponding e-folding times, 

2 £^puls MnO [Zn (a) ~ MiO (a)/'T'fco] ^0*^' 



2^2 
60 



Tsh = ---^ ' (54) 

(-Bpulssh) 11f^koY„kZk(a) 

^ _ 2i?puls W^^nO [2^;i(a) — MiO Kli Zi(a)/"6o] ^o«bo j-gg-j 

(^^puis b) Xnl^fiko Ynk Zk (a) (i9<^Ai/ dne)^ 



Now, using the asymptote ([29| for 1^^ and the asymptotes for uj and i?o, derived in Sec. IV, one gets for Tsh and tb 
in the limit T°° T^^(O), 

d2 / T^OO\" 

« ^^1- , (56) 

Up'R^ 

Tb cx finite. (57) 

Here a = —1 for the first model, a = 1/2 for the second, a = —1/2 for the third, and a = for the fourth model. The 
damping time Vo tends to some finite valu e for all the m odels. For the third model, presented in Fig. 7, the damping 
rate due to the shear viscosity drops as ^Jl — T°°/T^ at T°° — r^^(O) (that is, Tgh increases). Hence, very close to 
r^^(O), 1/Tsh <C l/rb so that r = rshTb/(Tsh + Tb) is bounded from above by Tb (this region of T°° is too narrow to 
be visible in the figure). 



VII. SUMMARY 



This work analyses the effects of finite temperatures on the oscillations of superfluid NSs. The most simple case 
of a radially pulsating nonrotating star is considered in the frame of general relativity. The impact of two temper- 
ature effects on the pulsations is examined. The first effect concerns the temperature dependence of the relativistic 
entrainment matrix Yi^ - one of the important ingredients of superfiuid hydrodynamics. The second effect consists 
in decreasing of the size of superfiuid region with increasing temperature. The first effect was studied in GA06, the 
second - in Refs. [H, H, HI] . 

We considered four different models which allow for one or both of these effects, and employ various profiles T^{r) 
of neutron critical temperature. We calculated pulsation spectra for these models. For each model the presence of two 
distinct classes of oscillation modes, the so called normal and superfluid modes, is revealed. The frequencies of normal 
modes almost coincide with the corresponding frequencies of a nonsuperfluid star. On the contrary, the frequencies 
of superfluid modes strongly depend on temperature and show different type of behavior depending on a model. It is 
demonstrated that the model which takes into account only the second temperature effect (dependence of the size of 
superfluid region on T) but ignores the flrst effect, yields qualitatively incorrect spectra. 

Using the approach developed in Refs. [H, [sj], we found the asymptotes for eigenfrequencies at T°° — > T^{0), 
which are in a good agreement with the results of numerical calculations. 

In addition, we analyzed the eigenfunctions for the first six oscillation modes. It is shown that the baryon current 
density is almost unperturbed in superfluid modes. Moreover, these modes are mainly localized in the core and 
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practically do not appear on the neutron-star surface. As a consequence, it might be difficult to detect them by 
observing modulation of the electromagnetic radiation from the star. 

We also calculated the characteristic damping times t for radial pulsations due to the shear viscosity and nonequi- 
librium modified Urea processes and analyzed the asymptotes of r at T°° — >■ T^{0). It is shown that damping times, 
calculated using the ordinary (nonsuperfluid) hydrodynamics may differ severalfold from t calculated for normal 
modes (even though the eigenfrequencies of normal modes almost coincide with the corresponding eigenfrequencies 
of a nonsuperfluid star). It is also demonstrated that, generally, damping of superfluid-type pulsations occurs by 1-3 
orders of magnitude faster than damping of normal-type pulsations. 
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